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Analytic methods for evaluating pointing errors caused by external disturbance 
torques are developed and applied to determine the effects of representative values 
of wind and friction torques. The expressions relating pointing errors to disturbance 
torques are shown to be strongly dependent upon the state estimator parameters, 
as well as upon the state feedback gain and the flow-versus-pressure characteristics 
of the hydraulic system. Under certain conditions , when control is derived from an 
un corrected estimate of Integra/ position error , the desired Type II servo properties 
are not realized and finite steady-state position errors result. Methods for reducing 
these errors to negligible proportions through the proper selection of control gain 
and estimator correction parameters are demonstrated. The steady-state error pro- 
duced by a disturbance torque is found to be directly proportional to the hydraulic 
internal leakage. This property can be exploited to provide a convenient method of 
determining system leakage from field measurements of estimator error, axis rate, 
and hydraulic differential pressure. 


I. Introduction 

Recent studies of mechanisms contributing to limit- 
cycle behavior have led to a need for more accurate mod- 
eling of the disturbance torque response characteristics of 
the 70-m antenna axis servos. Of particular interest in 
the limit-cycle studies is the transient behavior of the var- 
ious plant and estimator states that occur during friction- 
induced limit cycling. The traditional assumption that the 
estimator states accurately track those of the plant be- 
comes inaccurate when the plant is subjected to external 
disturbance torques. This shortcoming necessitated the 
development of the present multivariable axis servo model 
where the plant and estimator states are distinct. 

Precise pointing of the 70-m antenna is accomplished 
through the use of a two- axis, azimuth-elevation, bullgear- 


driven servo system. Control torques are produced by 
fixed-displacement, axial-piston hydraulic motors, which 
are coupled to the axis bullgears through spur-gear reduc- 
ers. Four such motor /gear reducers are employed for each 
axis. Backlash is eliminated by separate counter torque mo- 
tors which apply a constant torque bias to the output pin- 
ion of each control motor. The hydraulic connections to 
the countertorque motors are arranged as shown in Fig. 1 
so as to preload all four gear reducers and apply zero net 
torque to the bullgear. Torque modulation is accomplished 
by four port-hydraulic servo valves. Servo control con- 
sists of a hardware rate loop with tachometer feedback and 
a computer-based position servo employing state-variable 
feedback. 

The servo model is based on the nonlinear orifice-flow 
equation of the valve, along with the motor-flow equation 
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as described in [1]. A piecewise linear representation of 
the valve is obtained by partial differentiation of the valve 
equation (Eq. 5 of [1]) to yield the flow Q v as a linear 
function of the valve input current I v and load pressure 

Pl: 

Q v = K p I v -D h Pl (1) 

with 

K P = I< V P X J 2 


and 

P v = P s -Pl- Pr 

where Ps and Pr are the regulated system supply pres- 
sure and return pressure, respectively. As shown later, the 
load pressure is sufficiently small relative to the supply 
pressure such that the range of the valve pressure remains 
within a 2 to 1 ratio. This justifies the use of a constant 
value approximation for the flow gain Kp for control dy- 
namic analysis. In contrast, the equivalent damping Dr 
is linearly proportional to the flow, which varies in pro- 
portion to antenna rate. Therefore, Dr varies between a 
minimum value corresponding to leakage flow at the zero 
antenna rate, and a maximum corresponding to the maxi- 
mum tracking rate of the antenna. 

The hydromechanical system model incorporates the 
linearized valve of Eq. (1) along with a motor coupled to a 
rigid-body inertia load representation of the antenna struc- 
ture, as shown in block diagram form in Fig. 2, where Cr , 


Jm, and Vm represent the hydraulic compressibility, motor 
inertia, and motor displacement, respectively. 

The block diagram of the equivalent plant (servo-loop 
hardware) in Fig. 3 incorporates the model of Fig. 2, the 
tachometer feedback and control amplifier with associated 
compensation networks, and two additional integrations 
which produce the angular position and position integral 
states. The two inputs represent, the electrical rate com- 
mand input to the plant and an equivalent external dis- 
turbance torque; Kr represents a constant with value pro- 
portional to the rate loop gain, and P 2 and Z 2 are the pole 
and zero frequencies of the rate loop compensation net- 
work. Figure 3 also includes a simplified equivalent of the 
tachometer feedback network obtained by neglecting the 
network pole, which is at a relatively high frequency. This 
approximation introduces the acceleration feedback branch 
shown in Fig. 3 where the parameter Z\ corresponds to the 
negative real frequency of the network zero. 

To provide additional insight into the effects of rate 
and position loop parameters on system compliance, com- 
pliance equations are developed separately for the open- 
position loop case, for the hardware position proportional, 
integral, and derivative (PID) feedback, and for the closed- 
loop state variable controller. It is shown that the combi- 
nation of plant-state and estimator-state feedback in the 
precision mode leads to compliance characteristics different 
from those of the computer mode, even when both modes 
have identical plant and estimator dynamics. 

il. Compliance Equations for the Open 
Position Loop Case 

The open position loop compliance can be derived 
from the angular position response to a unit-step torque 
input. By application of the Mason transmittance rule to 
the block diagram of Fig. 3, the compliance transfer func- 
tion becomes: 


0m (1/Ja/s 2 )(1 + Dh /Chs){\ + Pi/s) 

Tx (1 4* P 2 / 5 )(l T £)// /Crs -f- Rr/CrJmZ.s + {Vjfa 4- K r) / J m C r s 2 ) 4- (Z 2 ~~ P 2 ){Rr! C r J m s 2 )(l / s 4- \jZ\) 


which leads to 


0 a/ (1/gXf 4- Dh/Cr)(s 4- P 2 ) 

Tx Jm(s 4* Pt)(s 2 4- ( Dr/Cr)s 4* Vm/JmCr ) 4- ( I<r/CrZi)(s 4- Z x )(s 4- Z 2 ) 
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Application of the final value theorem to Eq. (2) in- 
dicates a constant steady-state rate in response to a step 
function torque disturbance. The steady-state compliance 
is given by 

9m ^h_ 

T X ~ s(V* + K r Z 2 /P 2 ) 


III. Compliance Equations for Hardware FID 
Feedback 

The compliance of the closed-loop PID feedback sys- 
tem can be calculated by the same method used earlier for 
the open-loop case. The addition of position, integral, and 
rate feedback branches produces three additional denomi- 
nator terms in the compliance expressions. Thus, 


©Af (1 /Ja/s 2 )(1 4- Dh/Chs){1 T 

Tx (1 4- P 2 /s)(l 4- Dr/Crs 4- Kr/ ChJmZ i s + (Vm T Kr)/ JmCrs 2 ) 

4- (Z 2 - P 2 ){Kr/C h Jms 2 ){1/s + 1/Zi + Ki/s 3 + K 2 /s 2 + K 3 /s) 

4 - (1 4 - P 2 / s){K r/Ch J m s)(K\ / s 3 4 - K 2 /s 2 4 * K 3 /s) 


which leads to 


&m s(s 4- Dr/Ch)(s 4* P 2 ) 

17 ~ iZKs + a »)(* 2 + + v^/jmc h ) 

4 -(Kr/Cr)(s 4- Z 2 )(s 2 4- 5 3 jZ\ 4- 1\\ 4- K 2 s 4* A 3 S 2 ) 


( 3 ) 


and the steady-state compliance becomes 


Qm _ Drs 
T x ~ I<iKr(Z 2 /P 2 ) 


The compliance transfer functions of Eqs. (2) and (3) in- 
clude numerator zeros corresponding to the network pole 
P 2 and to the effective hydraulic damping Dr/Cr. The 
denominators are seen to contain the poles of the respec- 
tive open/closed-loop system. From the dependence of Dr 
on hydraulic flow, the steady-state compliance properties 
are seen to be antenna-rate dependent. For typical closed- 
position loop parameter values, the damping Dr has rel- 
atively little effect on the closed-loop poles. Therefore, 
since the shape of the compliance transient is determined 
by the locations of the transfer function zeros relative to 


the poles, the closed-loop transient properties are also seen 
to vary as a function of Dr . 

IV. Compliance Equations for the State Vari- 
able Controller 

Earlier methods of disturbance torque effects analy- 
sis [2] were based on the assumption of negligible errors 
in the estimator states in relation to the corresponding 
plant states. This permitted a simplification of the sys- 
tem model by substituting plant-state feedback in place of 
the estimator-state feedback, thereby reducing the number 
of states required in the model. This assumption, widely 
used in evaluating command input transient responses, was 
found to produce finite errors in the determination of dis- 
turbance transient responses of the axis servos. It was 
subsequently replaced by a superior method employing full 
modeling of the estimator states as well as the plant states. 
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x — Ax -f B U 


The model for disturbance torque response of the 
state- variable axis servo controller is based on a linearized 
multi-input state- variable representation of the system of 
Fig. 3, with the addition of the state-estimator and control- 
feedback gain. The plant state is represented by the gen- 
eralized state equations, where x and Y are the state and 
output by 


Y = Cx + DU 

respectively. 

The corresponding A, B, C, D matrices from Fig. 3 

are: 


A = 


0 1 
0 0 
0 0 
0 0 
0 0 


0 

0 

1 Jm 


—Vm/Ch — Kr/Ch —Dh/Ch — K r/ C ji J m Z\ Kji(Z 2 — P 2) /Cm 
-1 -IJJmZx -P 2 


B 


0 

0 

0 

K r /C h 

1 




D = 


0 

0 

1 /Jm 
— I\r/Ch 
-1/JmZx 


0 10 0 0 

'0 0 " 

0 0 
0 0 
0 0 


0 0 


Representing the estimator state by x , and estimator out- 
put by y, the estimator equations become: 

* = a* + B 1 U 1 + l(y-y) 

y = Cx^-DxUi 

where B 1 , B 2y Du D 2l U\ , and U 2 represent the first and 
second columns of B and D and the first and second ele- 
ments of U, respectively. This distinction is essential be- 
cause both inputs couple directly into the plant, while only 
the rate command input U\ couples directly into the esti- 
mator. 


V. Computer Control Mode 

Because the rate command U\ is formed differently for 
the precision mode, the two control modes are addressed 
separately. Incorporating the control-feedback gain K into 
the expression for the rate command input U\ for computer 
mode, and substituting into the plant and estimator equa- 
tions leads to 

Ui = —Kx 

x = Ax — B\Kx -f B 2 U 2 (4) 

x =: LCx + (A — BiK — LC — LD\K)x + LD 2 U 2 



With the estimator error, 


prising the plant and estimator states leads to the matrix 
differential equation form 


x — x — x 

x = (A - LC)x + LD\Kx + (B 2 - LD 2 )U 2 



A 

LC 


— B\K 

A -BiK-LC 



Ui (7) 


Since for the present case both D\ and D 2 are zero, the 
estimator and estimator-error equations simplify to 

x = LCx + (A - B X K - LC)x (5) 

and 

i = (A - LC)x + B 2 U 2 (6) 


which is compatible with existing linear system analysis 
and simulation software tools. The steady-state distur- 
bance torque compliance properties of the combined plant 
and estimator system can be determined from the steady- 
state solution of Eq. (7). The general expressions for the 
steady-state values of the individual plant and estimator 
states can be determined by a symbolic expansion of the 
determinants resulting from application of Simpson’s Rule 
to Eq. (7). 


The equations for the plant and estimator are illus- Since, in the steady state, the derivatives represented 

trated in the state-space block diagram in Fig. 4. Com- by the left-hand side of Eq. (7) equal zero, steady-state 

bining Eqs. (4) and (5) to form a single state vector com- position estimate x 2 becomes 



' 



1 1 

I I 




A | 

I 

-B\ki 

1 1 

1 b 2 1 

I I 


B\ &3 &4&5 

— det 

1 

1 

1 

— 

— 

1 i 

1 1 

1 o | 

— 



1 

A 

-1 

I 0 | 




LC | 

-Biki 

-LC 1 

1 0 1 

a 3> 

5 “ t 5 ] 


1 



1 0 | 







j 0 | 


_ 




A 

—BiK 





det 

LC 

A - BiK - 

LC 



( 8 ) 


where the numerator matrix is obtained from the 10 x 10 
denominator matrix by replacing the column correspond- 
ing to x 2 (column 7) with the right-hand-side vector, in 
accordance with Simpson’s Rule. The subscripts applied 
to capital-letter symbols designate the respective columns 
of the associated matrix, and A 3j 5 denotes a matrix com- 
prised of the third through fifth columns of A. 

A row by row examination of the numerator matrix re- 
veals that the sixth row is comprised of the five elements of 


the first row of TC, a single element (an — &n&i — and 

the last three elements of the first row of [A — B\K — LC], 
It will be seen that, since the first element of L is zero in the 
current parameter set, and because of the sparseness of A 
and B, the elements of the sixth row are all zeros and the 
numerator determinant vanishes. This indicates that, in 
the presence of a constant disturbance torque, the steady- 
state position estimate x 2 is identically equal to zero. This 
in turn indicates that the final values of the position x 2 and 
the position estimation error x 2 are identical, and x 2 can 
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thus be evaluated using Eq. (6). This approach avoids the 
complexity of evaluating the determinant of the 10 x 10 
matrices of Eq. (8). 

Thus, in the steady state, using Eq. (6), 
x = (A - LC)x + B 2 U 2 = 0 

and 


— det 

Ax 

b 2 a 3 A 4 A5 

det 

A -LC 



Note that L and C are absent from the numerator because 
all elements of C, except for the second, are zero. Note also 
that the first row of A and the first element of C are zero, 
causing both determinants to vanish. However, since the 
plant integral state is uncoupled from both the plant and 
the estimator, the corresponding first rows and columns 
can be deleted from both the numerator and denominator 
with no loss in generality. It should be noted that if the 
first element of C is assigned a nonzero value, the denom- 
inator determinant remains finite as long as L is nonzero. 
The implications of this integral position feedback will be 
discussed later. 

Substituting the values of A, Bn, and C into the ex- 
pression for Xn } 


0 0 

1/Jm 0 

—Dh/Ch — Kr/Ch Jm%i Kr(Z 2 ” Pt)/Ch 
-1/ JmZi —Po 

u 2 -1 2 10 0 

— / 3 0 1 /Jm 0 

det 

-/ 4 -VZJCh - I<r/C h -Dh!C h - Kr/C h JmZ\ K r {Z 2 — P 2 )/Ch 
-h -1 -1/ JmZi P’1 

Expanding the determinants leads to the nonzero steady-state position estimate error x 2} and since the steady-state 
plant position x 2 equals x 2 , 


0 1 

— 1 /Jm 0 

Kr/Ch JmZi -Vm/Ch - Kr/Ch 
1/JmZi -1 


Xn PfT 

U 2 ^2 iyjij + P rZ 2 / P 2 ) ~b ^3 (BjiJm T KrZ 2 /Z 1 P 2 ) + I^Ch + h(Z 2 /P 2 — 1) 


which, using current 70-m antenna servo parameter values, 
can be approximated by 

#2 ^ DliPl/Zn 

U 2 Kr \ 2 

This result invalidates the original premise that control 
feedback of estimated integral error is equivalent to feed- 
back of plant integral error in imparting Type II servo 
performance. That premise is valid only for inputs, such 
as the position command, which are coupled equally to the 
estimator and the plant. From Eqs. (3) and (4) and the 


block diagram of Fig. 4, it is seen that the plant is influ- 
enced by the control input U\ and a disturbance input {7 2j 
while the estimator is influenced by the same control input 
and by the estimator feedback error LCx . In equilibrium, 
the plant disturbance U 2 is compensated by the control in- 
put U\\ thus, the identical U\ input to the estimator must 
be counteracted by the error LCx in order to obtain es- 
timator equilibrium. This implies that a finite estimator 
error x will always result from a plant disturbance input, 
while the individual components of x will be determined 
by the product of L and <7. This line of reasoning explains 
the absence of the control gain K from the steady-state 
compliance expression. 
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In the present case, where there is no feedback of 
x\ } the integral estimate error Xj is allowed to grow un- 
bounded. The addition of a small amount of plant integral 
feedback to the estimator by assigning a small value to c\ 
would bound the error xi, thereby forcing X 2 , £ 2 , an d £2 
to zero. The required nonzero x would then result from 
nonzero values of the components other than and X3. 
This additional feedback imparts the desired disturbance 
accommodating control [3] properties to the system. 


VI. Precision-Control Mode 

In the precision (autocollimator feedback) mode, the 
control input U 2 is derived from the autocollimator, an 
electro-optical device that senses plant position error di- 
rectly. This error signal is filtered to remove noise and 
then integrated. The resulting plant integral and position 
errors are then combined with additional damping terms 
provided by the position estimate to form the plant input 
according to 


the steady-state estimator error x is therefore identical in 
both the computer mode and precision mode. However, 
the altered form of the control gains Kp and Ke in the 
precision mode change the coupling of the integral position 
xi and its estimate Xi. Accordingly, when solving Eq. (9) 
for final values, the row and column corresponding to the 
integral estimate X\ are deleted, and those corresponding 
to the integral x\ are retained. As a result, the final value 
solution of Eq. (9) yields the expected zero value for the 
position X 2 and a nonzero position estimate x 2 . This result 
implies a nonzero estimator error x consistent with Eq. (6). 

These observations demonstrate an interesting duality 
between the two control modes, whereby the substitution 
of plant states for their corresponding estimator states in 
the control gain law results in the interchange of the fi- 
nal values of plant and estimated position variables. This 
tends to confirm the validity of the general expressions of 
Eqs. (7) and (9), and indicates that in the precision-control 
mode, a zero steady-state error always results from a dis- 
turbance torque input. 


U\ — — £ 1 X 1 — £ 2 X 2 — — k 3 x 4 — fc 5 x 5 


The resulting system equations therefore take a slightly 
different form from that of Eq. (7) due to the mixed plant 
and estimator feedback. Thus, for precision mode, 


x 

x 


A-BiKp -B x K e 
LC - B\K P A - BJ< e - LC 


x 

X 


VII. Compliance Transient Properties 

The shape of the transient response to disturbance 
torques can be inferred from the relative locations of the 
poles and zeros of the compliance transfer functions for 
the various system configurations. For the open-loop and 
hardware PID feedback systems, the poles and zeros are 
defined by Eqs. (2) and (3), respectively, where it is seen 
that the poles are the same as those associated with the 
respective command input-output transfer functions. 


+ 


B 2 — BiKpD 2 

0 


u 2 


( 9 ) 


where Kp and Kp are the gain vectors associated with the 
plant and estimator states respectively. For the precision 
mode, 


For the state variable controller configurations, the 
poles are the respective eigenvalues of the square system 
matrices in Eqs. (7) and (9). The corresponding zeros are 
the complex frequencies of zero response of x, x, and Y for 
any disturbance input U 2 and initial condition. Thus, for 
the computer-mode configuration of Eq. (7), the response 
zeros are the roots of 


' K P ' 


k\ k 2 0 0 0 

K e 

L J 


0 0 k 3 k A k 5 


The B\KpD 2 term is included in Eq. (9) to maintain gen- 
erality, even though the KpD 2 product is always zero in 
the absence of plant acceleration feedback. 

From Eq. (9), it is seen that Eq. (6) for the estimator 
error is also applicable to the precision mode case, and 


(si — A) —B\K —Bo 


X 


‘O' 

LC (sI-A + BJ< + LC) 0 


X 

— 

0 

) 

O 

O 

O 


U 2 _ 


_ 0 _ 


To avoid the longhand expansion of the above determinant, 
the zeros were determined by a numerical evaluation of 
the transfer function zeros of the system of Eq. (7) with 
coefficient values representative of the 70-m azimuth servo. 
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The results coincide precisely with the two nonzero real 
roots of the determinant 

[si — A] B 2 
C 0 

at - Dh /Ch and — P 2 (the root at zero is absent), and with 
the five eigenvalues of [A - B\ K - LC ]. A numerical com- 
putation of the eigenvalues of the square matrix of Eq. (7) 
shows that the accompanying poles are the five eigenvalues 
of [A — BiK ] and the five eigenvalues of [A - LC]. 

For the precision-mode configuration of Eq. (9), the 
response zeros are the zeros of 

[si — A + B\ Kp] — B\I\e —B 2 

LC-BxKp [%1 -A + B x K e + LC\ 0 
C 0 0 


that the initial transient characteristics of the state vari- 
able controllers are governed solely by the plant parame- 
ters, and the departure from this initial characteristic is 
governed by the controller. 

VIII. Numerical Results 

Numerical evaluations of the differential equations in- 
corporated the physical parameter values of Table 1, the 
control coefficients of Table 2, and the disturbance torque 
parameters of Table 3. Physical parameter values of Ta- 
ble 1 were derived from component specifications accord- 
ing to [1] with the rate loop gain and network parameters 
based on [4]. The Kr value of 5.107 corresponds to a rate 
loop dc gain product of 40. The values shown for damp- 
ing Dr are calculated using Eq. (1) and correspond to the 
axis rate range of 1 to 40 times the sidereal rate with added 
leakage Qhl. Aerodynamic parameters used in determin- 
ing wind torques were extrapolated by McGinness 1 from 
earlier wind-tunnel tests performed on scale models of the 
64-m antenna. 



X 


'o' 

X 

X 

— 

0 




_o 


Numerical evaluation of the transfer function zeros indi- 
cates two real zeros at — Dr/Ch and — P 2 , two zeros 
at the origin (one of which is canceled by a pole), and 
four additional real and complex zeros which are near, 
but noncoincident with, four of the poles of the square 
matrix of Eq. (9). This result indicates that unlike the 
computer-mode case of Eq. (7), the precision mode ex- 
hibits zero steady-state disturbance torque compliance. In 
addition, the low-frequency zero from the eigenvalues of 
[A — B\K — LC] arising from Eq. (7) is absent in this 
case, resulting in a faster disturbance recovery transient 
in the precision mode. The system poles computed from 
Eq. (9) are identical with those from Eq. (7) for the com- 
puter mode. 

From the foregoing and from Eqs. (2) and (3), it is 
seen that the real zeros at —D h /Ch and -P 2 appear in 
the compliance expressions for the open-loop and PID feed- 
back systems, as well as for the computer and precision- 
mode state variable controllers. The differing steady-state 
compliance properties of these four configurations arise 
from the presence of poles or zeros at the origin. This 
commonality implies that all four configurations exhibit 
identical disturbance transient characteristics for a short 
time interval following the transient. This in turn implies 


The position feedback control gain K and estimator 
gain L shown in Table 2 were calculated according to the 
optimal control criteria of Alvarez and Nickerson [5]. The 
baseline values L\ and K 4 in the table were calculated to 
duplicate the dominant closed-loop poles of [5] when ap- 
plied to the model of Fig. 3. Alternate parameter sets 
were calculated in a similar manner with weighting ad- 
justed to shift the dominant poles inward or outward from 
the origin. Those estimator gains that include integral 
terms were calculated using an alternate output vector 
H = [1 1 0 0 0], which included equal weighting of 

integral and position. 


The wind torque moments in Table 3 were calculated 
according to McGinness 2 using the torque relationship 


torque = 


CrgrP^ 

4 


(10a) 


where Ct is the aerodynamic lift coefficient, the dynamic 
wind pressure is 


<7 = 



(10b) 


with air density, p — 0.00238 lb/ft 3 , and wind velocity, v. 


1 II. McGinness, “Effects of Wind Loading on 64- And 72- 
Meter Diameter Antenna” (JPL internal report), Reorder 
No. 84-2, May 1984. 

2 McGinness, ibid. 
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From Eqs. (10a) and (10b), the disturbance torque is 
seen to be a quadratic function of the wind velocity, so the 
disturbance imparted by a step change in velocity depends 
on the initial and final velocities rather than the step am- 
plitude. Thus, if v m is the mean of the initial and final 
velocities and Vd is their difference, the torque disturbance 
is proportional to the product v m Vd. The worst-case dis- 
turbance then results when a large step- velocity change is 
superposed on a large mean velocity such that their prod- 
uct is a maximum. Assuming 30 mph as the maximum 
average wind and 12 mph as the maximum gust, the worst- 
case disturbance corresponds to a step change between 24 
and 36 mph. 

Breakaway friction torque levels were extracted from 
acceptance-test data recorded when the overseas antennas 
were built in 1972. The range of values shown is consistent 
with more recent informal reports of observations at the 
Goldstone site. 

The differential hydraulic pressures corresponding to 
the disturbance torques are included in Table 3. Because 
of its linear relationship to torque and its ease of direct 
measurement, the differential pressure has become a fa- 
miliar unit of torque measurement to those working with 
the antenna. 

The disturbance torque transient responses of the 
open-loop, computer-mode, and precision-mode configu- 
rations obtained from the time solutions of Eqs. (7) and 
(9) with I\ = I< 4 and L — L\ are shown in Fig. 5. The 
disturbance input corresponds to a 12-mph wind step com- 
bined with a 30- mph mean wind. The open-loop case was 
obtained from Eq. (7) with zero value of control gain K, 
The comparatively rapid recovery of the precision mode 
and the finite steady-state error of the computer mode are 
clearly visible in the time responses. The slower recov- 
ery of the computer mode is attributed to a low-frequency 
transfer function zero arising from Eq. (7), which is absent 
from the precision-mode case. 

Figure 6 shows the computer-mode disturbance torque 
transient response for various antenna rates. The re- 
sponses were generated from Eq. (7), where the Dh/Cjj 
term in the A matrix was adjusted according to Eq. (1) for 
each antenna rate. The effect of the changing Dh /Cjj zero 
in the transfer function is evident from the change of the 
transient shape and of the final value. Figure 6 can also 
be used to judge the effects of increased hydraulic leakage 
by converting a known volumetric leakage to an equivalent 
antenna rate through the hydraulic displacement and gear 
ratio. The characteristic corresponding to the sum of the 


leakage equivalent rate and actual rate then represents the 
system behavior. 

Figure 7 shows the computer-mode transient response 
for the various values of control gain K listed in Table 2. As 
indicated earlier from the general properties of Eqs. (6) and 
(7), the gain K is seen to influence the speed of recovery, 
but has no effect on the final value. 

Figure 8 shows the computer-mode transient response 
for the various values of estimator gain L listed in Table 2. 
The effect of those L vectors having finite integral error 
correction terms l\ is seen in the faster recovery time, as 
well as reduction of the final value. The final value for 
those cases is still finite due to the absence of position 
integral feedback. 

Figure 9 shows the precision-mode transient response 
for the various values of control gain K listed in Table 2. 
Increasing values of K are seen to decrease the peak error 
and speed the recovery transient. 

Figure 10 shows the precision-mode transient response 
for the various values of estimator gain L listed in Table 2. 
As expected, the estimator gain L has little effect on the 
compliance transient when the response of the estimator 
is faster than that of the control gain. 

IX. Summary and Conclusions 

The transient disturbance-torque rejection properties 
of the 70-m axis servos are shown to be governed by phys- 
ical hardware parameters, as well as by the properties of 
the software-control algorithm. In particular, the effective 
hydraulic damping Dh/Ch> which is strongly affected by 
hydraulic leakage as well as by axis rate, determines the 
peak transient error, and in the case of the computer mode 
it also produces a finite steady-state error. The peak and 
steady-state disturbance-torque errors are essentially un- 
affected by the control gain K } but can be reduced by 
increased estimator gain L at the expense of noise rejec- 
tion. 

The precision-control mode, by virtue of the direct 
feedback of hardware (as opposed to estimated) position 
and integral errors, is shown to have disturbance torque re- 
jection properties superior to those of the computer mode. 
This deficiency of the computer mode results from the 
omission of the position integral error from the estimator 
equations, and can be corrected by modifying the equa- 
tions to include the integral. Significant improvement can 
also be obtained by modifying the estimator gain parame- 
ter L to include a coefficient corresponding to the integral 
position estimate. 
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Table 1. The 70-m antenna axis servo parameter values 


Symbol 

Pair ame ter 

Units 

Value 

Jma 

Inertia moment , azimuth 

in.-lb-sec 

12.0 

J ME 

Inertia moment, elevation 

in.-Ib-sec 

8.0 

V H 

Hydraulic displacement 

in. 3 /radian 

1.53 

Ch 

Hydraulic compressibility 

in. 3 /psi 

0.00314 

D h 

Hydraulic damping 

in. 3 /psi-sec 

0.00188 to 0.0754 

Qhl 

Hydraulic leakage 

in 3 /sec 

1.67 

N 

Reduction gear ratio 

dimensionless 

28730 

Pi 

Tach network negative pole frequency 

sec” 1 

80 

P 2 

Lag network negative pole frequency 

sec” 1 

0.24 

Zi 

Tach network negative zero frequency 

sec” 1 

5.0 

z 2 

Lag network negative zero frequency 

sec” 1 

4.4 

Kr 

Rate loop gain constant 

in. 6 

5.11 

K V 

Valve flow constant 

in 3 /mA-sec/psi°- 5 

0.731 

Pv 

Valve pressure drop 

psi 

2750 

Cta 

Aerodynamic lift coefficient, azimuth 

dimensionless 

0.1217 

Cte 

Aerodynamic lift coefficient, elevation 

dimensionless 

0.1231 

D a 

Reflector aero dynamic diameter 

m 

70 



Table 2. System, control gain, and estimator parameters 


(a) System matrices 

A = 

B = 

C = 
D = 

(b) Control gains K 
Control gain K 

I<i 

Ki 
Kz 
K 4 
K b 


0 

0 

0 

13.2484 

—0.2400 


0 0 

0 0 

0 0.0833 

16.2643 -0.2711 

5.1070 -0.0851 


0 1 .0000 0 0 0 


0 1.0000 
0 0 

0 0 

0 0 

0 0 


0 

1.0000 

0 

-23.7194 
- 5.1070 


0 

0 

8.3333 
-27.7072 
- 8.5117 


0 


Numerical value 


Control poles from eigenvalue (A — BK) 


0.5200 

1.3417 

0.5712 

0 

0 

0.3162 

0.9097 

0.1685 

0.0162 

0.1416 

0.4472 

1.1113 

0.2087 

0.0279 

0.1716 

0.6325 

1.3684 

0.2642 

0.0500 

0.2106 

0.8944 

1.7032 

0.3450 

0.0916 

0.2629 

1.2649 

2,1515 

0.4690 

0.1699 

0,3339 


-12.2372 4- 12.5297* 
-12.2372 - 12.5297* 
-2.5455 

-0.4636+0.4268* 
-0.4636 - 0.4268* 

-12.8664 + 9.5595* 
-12.8664 - 9.5595* 
-2.4141 

-0.3937+0.3861* 
-0.3937- 0.3861* 

-12.9520 + 9.4886* 
-12.9520- 9.4886* 
-2.4404 

-0.4668 + 0.4539* 
-0.4668 - 0.4539* 

-13.1214 + 9.3446* 
-13.1214 - 9.3446* 
-2.4911 

-0.5509 + 0.5291* 
-0.5509 - 0.5291* 

-13.4526 + 9.0474* 
-13.4526- 9.0474* 
-0.6445 + 0.6078* 
-0.6445 - 0.6078* 
-2.5861 

-14.0877+ 8.4129*' 
-14.0877- 8.4129* 
-0.7426 + 0.6822* 
-0.7426 - 0.6822* 
-2.7556 


Table 2. System, control gain, and estimator parameters (contd) 


(c) Estimator gains 
Estimator gain L 

Numerical value 

Estimator poles from 
eigenvalue (A — BK\ — LC) 

Compliance zeros from 
eigenvalue (A — BKa — LC) 

L x 

0 

-12.7800 + 9.6298* 

-13.1060 + 9.3382* 


0.4546 

-12.7800 - 9.6298* 

-13.1060 - 9.3382* 


0.0033 

-2.3859 

-2.1842 


-0.0020 

-0.4559 

-1.4267 


0.0026 

0 

-0.4674 

l 2 

0.8638 

-12.4988 + 9.7652 t 

-12.6999 + 9.4375* 


3.1261 

— 12.4988 — 9.765 2t 

-12.6999- 9.4375* 


4.8336 

-4.8845 

-5.0267 


— 7.6308 

-1.1912 

-2.5190 


-9.1392 

0 

-0.0162 

Lz 

0.8075 

— 12.7656 + 9.6345* 

-13.0283 + 9.3205* 


1.4625 

-12.7656-9.6345* 

-13.0283 - 9.3205*' 


1.1140 

-1.9392 + 0.7044* 

-2.6041 + 1.3095* 


-0.6472 

-1.9392 - 0.7044* 

-2.6041 - 1.3095*' 


0.5716 

0 

-0.0333 

L a 

0.6204 

-12.7686 + 9.6345* 

-13.0829 + 9.3379* 


0.5545 

-12.7686 - 9.6345* 

-13.0829 - 9.3379* 


0.1357 

-2.5083 

-2.0488 + 0.4207* 


-0.3804 

-0.4562 

-2.0488- 0.4207* 


-0.5582 

0 

-0.1267 


Table 3. Typical disturbance torque levels 


Azimuth 

Elevation 

Wind moments for 30-mph wind 

2.67 x 10 6 ft-lbs 

2.70 x 10 6 ft-lbs 


729 psi 

738 psi 

Breakaway friction torque 

1.03- 1.28 X 10 6 

0.84 - 1.46 X 10 6 ft-lbs 


280 — 350 psi 

230 — 400 psi 
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Fig. 1. The 70-m antenna axis servo countertorque system. 



Fig. 2. Servo hydromechanlcai system linearized model. 
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Fig. 4. Servo position controller block diagram. 



Fig. 5. Disturbance torque transient response of open-loop, computer-mode, and 
precision-mode configurations. 
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Fig. 6. Disturbance torque transient response for various axis rates. 



TIME, sec 


Fig. 7. Computer-mode disturbance torque transient response for various K. 




AZIMUTH ERROR, mdeg AZIMUTH ERROR, mdeg 



Fig. 8. Computer-mode disturbance torque transient response for various L 
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TIME, sec 

Fig. 10. Precision-mode disturbance torque transient response for various L 




